arXiv:1501.01850v2 [gr-qc] 13 Jan 2015 


Ricci Time in the LemaTtre-Tolman Model and the Block 

Universe 


Yasser Elmahalawyf Charles Hellabyt and George F.R. Ellis^ 

Department of Mathematics and Applied Mathematics, 

University of Cape Town, Rondebosch, 7701, South Africa 


Abstract 

It is common to think of our universe according to the “block universe” concept, which says 
that spacetime consists of many “stacked” 3-surfaces, labelled by some kind of proper time, r. 
Standard ideas do not distinguish past and future, but Ellis’ “evolving block universe” tries to 
make a fundamental distinction. One proposal for this proper time is the proper time measured 
along the timelike Ricci eigenlines, starting from the big bang. This work investigates the shape 
of the “Ricci time" surfaces relative to the the null surfaces. We use the LemaTtre-Tolman metric 
as our inhomogeneous spacetime model, and we find the necessary and sufficient conditions for 
these {r= constant} surfaces, S(t), to be spacelike or timelike. Furthermore, we look at the 
effect of strong gravity domains by determining the location of timelike S regions relative to 
apparent horizons. We find that constant Ricci time surfaces are always spacelike near the big 
bang, while at late times (near the crunch or the extreme far future), they are only timelike 
under special circumstances. At intermediate times, timelike S regions are common unless the 
variation of the bang time is restricted. The regions where these surfaces become timelike are 
often adjacent to apparent horizons, but always outside them, and in particular timelike S regions 
do not occur inside black holes. 


1 Introduction 

There are numerous theories that discuss the nature of time and the way to measure time. Rovelli [37] 
and Barbour [3] assumed that time does not exist and does not “roll on”, according to the idea of 
standard block universe. Davies [11] proposed that time is real but flow is not. Price [36] and Barbour 
[3] argue that a block universe representation denies that time passes and embodies the idea that the 
flow of time is an illusion. Ellis [14] examined the various views and argued that time exists and its 
flow is real, calling this the evolving block universe (EBU). This idea allows us to contemplate the dif¬ 
ferent natures of the past and future, and differentiate between them in our spacetime representation. 
Moreover, it shows how to indicate the present “now” in the spacetime diagram. 


This paper will not repeat these discussions: the main issue we address is the concept proposed in 
[14] of how to characterize preferred time surfaces in a curved spacetime filled with a non-zero energy 
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density (for example the Cosmic Microwave Background Radiation, or CMB). This proposal can be 
investigated independently of its relation to the EBU idea, as it is an interesting geometrical question 
in its own right as to how the preferred spatial sections mentioned in [14] evolve. 

According to the EBU, the universe consists of a sequence of curved surfaces, labelled by proper 
time r, which is measured along a family of chosen timelines x l (v ) from the start of the universe to 
the present day. Ellis [14] argued that the present time surface at proper time r = r 0 is the surface of 
{r — r 0 } defined in this way, which can be measured by integration along a family of fundamental world 
lines; as time passes, r 0 continually grows, so the future boundary of spacetime is always extending 
(hence the name EBU). The fundamental world lines used in this definition are chosen to be the 
Ricci eigenlines, which represent the average motion of matter, and will be well defined in all realistic 
spacetimes, because of the existence of an all-pervading cosmic blackbody radiation in the universe. 

This prescription defines a family of surfaces that may be locally either spacelike, null, or timelike. 
The main idea of this paper is to examine the behavior of the {r = constant} surfaces S(r c ) and find 
out their nature: are they locally spacelike, timelike, or null? This is an interesting geometric question 
that can be asked independently of whether one accepts the idea of an evolving block universe or not. 
One might expect that these surfaces will be spacelike except in the case of very intense gravitational 
fields. We will show, by investigating these surfaces in a number of LemaTtre-Tolman models, that this 
is not necessarily the case. 

The LemaTtre-Tolman (LT) metric is an exact solution of the Einstein field equations (EFEs), 
with a dust equation of state, suitable for studying the possibility of an inhomgeneous universe, or 
inhomogeneous structures within it. Like the Friedmann-LemaTtre-Robertson-Walker (FLRW) model, 
it is spherically symmetric with comoving matter, but LemaTtre-Tolman models are generally inhomo¬ 
geneous in the radial direction. It was first proposed by LemaTtre [31] in 1933 and then by Tolman [44] 
in 1934; a few years later, in 1947, Bondi [6] investigated this model again. It has been widely used to 
model large and small scale structures, and to study gravitational collapse, cosmological observations, 
and exact nonlinear features of the EFEs [21, 23, 4, 2, 39]. Interestingly, it can be used to describe 
a homogeneous cosmology in one limit, and non-vacuum black holes in another limit [18, 29]. See 
[27, 35] for a summary of work on the LemaTtre-Tolman model, and [22] for a well illustrated outline. 


2 The LemaTtre-Tolman Model and its Evolution 

2.1 LemaTtre-Tolman Model 

The LemaTtre-Tolman (LT) metric is written in spherically symmetric, synchronous coordinates, 

ds 2 = -d t 2 + dr 2 + R 2 (d6 2 + sin 2 9 d</> 2 ), (1) 

where R is the areal radius, f(r) is a free function determining the local geometry, and R' = dR/dr. 
The matter is a pressure-free perfect fluid [16], comoving with the coordinates, and the energy mo¬ 
mentum tensor is defined by 


T ab = pu a u b , 

where p is the mass-energy density, and u a is the fluid’s four velocity, 

u a = 51 


( 2 ) 

( 3 ) 
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The solution of the EFEs gives the evolution equation 


R 2 


2 M{r) 
R 


f{r) + 


A R 2 


( 4 ) 


where R = dR/dt, the second free function M(r) is the gravitational mass within the comoving shell 
of radius r, and f(r) has an additional interpretation as / = 2 E(r) which is twice the local energy per 
unit mass of the dust particles. The EFEs also give the density as 


up 


2 M' 

R 2 R r 


The Kretschmann scalar [22] is an invariant measure of spacetime curvature which is given by 


( 5 ) 
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/C only diverges where R or R' are zero, while M and M' are not. Similarly, p only diverges where 
M'/R' is zero. The solutions of (4), when the cosmological constant A is zero, are of 3 types, depending 
on the value of / (or more correctly, the value of //M 3 / 2 ), 


Hyperbolic, / > 0 : 
Parabolic, / = 0 : 
Elliptic, / < 0 : 


R 

R 

R 


y- (cosh r/ - 1), 


M I — 


M 

FT) 


(1 — cos rj), 


(sinh r) — rj) 



(77 — sin 77) 


f 3 / 2 (t — a) ' 

M 
('t — a) 

M ’ 

(- f) 3/ \t- a l 
M 


( 7 ) 

( 8 ) 
( 9 ) 


where the last free function a(r) is the local time of the big bang, i.e. the time on each worldline at 
which R = 0. The bang time does not have to be the same for each worldline; usually outer shells of 
space and matter “explode” off the initial singularity at earlier times. 

The parameter 77 can be thought of as giving the stage of evolution, such as “early", “interme¬ 
diate”, and “late”. The 3 types of solution are not mutually exclusive, and it is possible to have a 
hyperbolic region outside and elliptic region, with a parabolic locus at the boundary between them; this 
models gravitational collapse in an expanding universe [28, 30]. Though extended parabolic regions are 
possible, they are all asymptotically homogeneous, consequently they will not be investigated explicitly 
here. The 3 arbitrary functions of this model allow a freedom to rescale the r coordinate plus two phys¬ 
ical freedoms, such as a(M ) and f(M). However, as noted in [22], each function contributes to the 
physical inhomogeneity, and specific “coordinate choices” restrict physical possibilities. For example, 
if M = r 3 is specified, then vacuum regions M' = 0 are excluded; similarly for the other two. 

A scale length and time, characteristic for each worldline, are defined as 

M ~ M 

R(r) = j Yy and T ( r ) = jm 72 ’ ( 10 ) 


respectively. Note that the lifetime from the bang to crunch in elliptic models is 2ttT. 


2.2 Singularities 

Apart from the big bang and big crunch singularities, which are non-simultaneous versions of their FLRW 
counterparts, LT models may also suffer from shell crossing singularities, as well as “shell focussing” 
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singularities. The big bang occurs on each worldline where t = a, causing R(a,r ) = 0 for all r. The 
big crunch occurs in elliptic models where t — a + 2-kT, also causing R = 0 on each worldline. The 
bang and crunch 3-surfaces are space-like [17] as will be seen in section 3, and have diverging density 
(5) and curvature (6). 

Shell crossings occur due to the collision between an inner and an adjacent outer shell of constant 
r, meaning R' — 0 locally; consequently density and curvature diverge. The shell-crossing surfaces 
are time-like, and they have a different redshift structure from the big bang[24]. Though matter 
collision may be allowable in a much more sophisticated matter description, shell crossings are physically 
unacceptable in LT models. The conditions on the LT arbitrary function that ensure no shell crossings 
are avoided were given in [24, 29]. See table 1. 

Shell focusing singularities can occur on the central worldline (origin), at the moment of the big 
crunch. Under certain conditions, many light rays can be emitted from this one point on the crunch 
and can even travel to infinity, rather than being captured by the crunch [26, 45, 34, 12, 9]. It will not 
be considered here. 


R 

/ 

M\ /', a 

> 0 

> 0 

M' > 0 
/ > 0 
a < 0 

but not all 3 equalities at once 

< 0 

M' > 0 

2tiT' + a > 0 
a < 0 

but not all 3 equalities at once 

= 0 

= 0 

o 

II ° ° 

"'+-T e 

< 0 

> 0 

M' < 0 
/ < 0 
a > 0 

but not all 3 equalities at once 

< 0 

M' < 0 

27 r T' + a < 0 
a > 0 

but not all 3 equalities at once 


Table 1: Conditions for no shell crossings. 


2.3 Origins and Regular Extrema 

An origin occurs where R(r 0: t ) = 0 for all f on a particular worldline, r 0 [7]. This implies that all 
time derivatives of R along r 0 are also zero. Obviously, even at the origin, t must vary over a finite or 
infinite range, and so must the parameter rj. So it is clear from equations (7) and (9) that, M/|/| 3 / 2 
must remain finite and non-zero in the limit as r —> r 0 , and M/\f \ must go to zero to make R zero. 
Therefore, M ~ |/| 3//2 in the neighbourhood of the origin, and both M and / go to zero there. But this 
does not mean that the time evolution is parabolic [17], When doing numerical calculations, we need to 
avoid calculating zero over zero, so for each choice of arbitrary functions, we need separate expressions 
for quantities like M// 3 / 2 that have well behaved origin limits. The combinations of arbitrary functions 
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and derivatives in equations (14)-(10) all require this treatment. In contrast to the big bang, an origin 
is timelike and non-singular. 

In a spatially closed spherical model, such as the k — +1 FLRW model, there must be place where 
the areal radius is maximum, decreasing on either side towards an origin; a closed spherically symmetric 
model has two origins. Similarly, in LT models of SKS black holes, vacuum and non-vaccum, there 
is a minimum of the areal radius at the “throat” or “neck”. Regular extrema in LT models occur if 
R' = 0 at one or more particular r m values without the density or curvature diverging and without a 
shell crossing forming [24], This means that 

R'(t,r m ) = 0 Vt (11) 

The conditions for a regular maximum or minimum [24] without shell crossings or surface layers are 

M'(r m ) = f'(r m ) = a'(r m ) = 0, f(r m ) = -1. (12) 

Unlike shell crossings, regular extrema remain at constant r, and are not singular. 


2.4 Special Cases 

The LT model has several interesting special cases; the first two are of particular interest here. The 
dust FLRW model is obtained when 


/ oc M 2/3 , a = 0. (13) 

If M' = 0, then we have at least part of the vacuum spherical spacetime. To get the full Schwarzschild- 
Kruskal-Szekeres (SKS) spacetime, requires an elliptic region, and a worldline where / = —1 and 
a' — f — 0 at which R' = 0, a' is maximum, and f and R' are minimum. One may easily construct a 
non-vacuum black hole with the same geometry and topology, if M' is made non-zero and minimum at 
/ = —1 [20, 18]. The Datt-Kantowski-Sachs metric is a well-behaved limit of the LT model, occuring 
if R — 0 and / = —1 globally [20]. The Vaidya metric is the limit in which the matter worldlines are 
null; after appropriate transformations, the limit of infinite local energy / —> oo is taken [32, 19]. 


2.5 Parametric evolution of R' 


Parametric expressions for the evolution of R' [24, 22] will be very useful for calculating the radial null 
slopes for both elliptic and hyperbolic cases in section 3. They are derived from (7) and (9): 


Elliptic, / < 0 : 
where 


Hyperbolic, / > 0 : 


R' 


4>i 

^2 

^3 


R' 


Tiyy^ 1 + “T^ 2 + (- a ')v / (-/)^3, 
2 + '7 sin// 

(1 — cos rj) ’ 

^ sin 77(377 — sin 77) 

2(1 —cos 77) 
sin 77 
1 — cos 77 


(14) 

(15) 

(16) 

( 17 ) 

(18) 
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where 


(19) 


^4 

^5 

A 


rj sinh rj ^ 

(cosh rj — 1) 

sinh rj (sinh rj — 3 rj) 
2 (cosh r/— 1) 
sinh rj 
cosh rj — 1 


( 20 ) 

( 21 ) 


The limiting behaviour of the evolution functions at early and late times, and their values at 
maximum expansion, will also be useful: 


Elliptic 

Hyperbolic 


Early 

Max Expansion 

Late 


Early 

Late 

V 

-)• 0 

= 7r 

5 = 27t — rj — > 0 

V 

0 

— > 00 

A 

~ ?7 2 /6 

= 2 

~ 47t/<5 

i>A 

~ rj 2 /6 

~ rj 


~ ?7 4 /40 

= 2 

~ 67 t/5 

^5 

~ 77 4 /40 

~ e ,? /4 

i >3 

~ 2/rj 

= 0 

~ —2/S 

A 

~ 2/rj 

~ 1 + 2e _,? 


Table 2: Limiting behaviours of the evolution functions. 


Finally, the rj derivatives of the hyperbolic ?/y functions are: 

dip 4 sinh rj — rj 

drj cosh /y — 1 ’ 

d ^5 sinh rj 3 (sinh rj — rj) 

drj 2 2(coshr/ — 1) ’ 

d^ 6 _ -1 

drj cosh rj — 1 

2.6 Scaled Conformal Time 

Below we will numerically calculate and plot results for a number of explicit models, and they will 
display a specific range of t and r values. The evolution calculations will be based on the conformal 
time rj, so it will be necessary to set a suitable range of rj and to have sufficient data points in that 
range to make a good plot. The relation between t and rj depends on the scale time T, via (7), (9) and 
(10). In addition, near parabolic regions, where / passes through zero, it can be seen that (7) and (9) 
imply rj —>• 0. To get a useful parameter in the parabolic case, the defintion rj = r]/y/\f \ is used before 
taking the limit / —» 0. These two problems are handled numerically by assuming a regular origin and 
defining the rescaled conformal time to be 


( 22 ) 

(23) 

(24) 


f 1/2 

«= 517 n - < 25 > 

3 Ricci Time Surfaces in LT Models 

In the context of the LT model, we now investigate under what circumstances the surfaces of constant 
Ricci time, S(t) are spacelike, null or timelike, for both hyperbolic and elliptic cases. This involves 
comparing the slopes of the S(t) surfaces with the slopes of the null rays. Furthermore, we derive the 
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necessary and sufficient conditions on the three arbitrary functions of the LT model, for the S(t) to 
be timelike at early, late and intermediate times. The results indicate the character of r and provide a 
guide to construction numerical examples. 


3.1 Ricci Time 


Ellis [14] argued that the “present” is a surface of constant time S'(r), and that the time r for any 
S(t) can be evaluated by integrating proper time along suitable worldlines, starting at the beginning 
of the universe up to S(t). He proposed that the paths to use are the Ricci eigenlines, that is, the 
integral paths x a {v ) of the timelike eigenvectors u a = dx a /dv of the Ricci tensor: 1 

T ab u b = Ai u a <=> R ab u b = A 2 u a , (26) 

pv 

x a (v)= / u a (v)dv (27) 

Jo 

pv 

t = / \/ \u a g ab u b \ dv. (28) 

Jo 

The resulting time r will be called the Ricci time. In the LT model, the Ricci eigenlines turn out to 
be the matter worldlines u a , as is evident from (2), and the proper time is just the local cosmic time 
since the bang, 


t = t — a. 

Consequently, the slope of S(t) is simply the derivative of (29), holding r constant, i.e. 

dt I , 



This already shows that the slope of the big bang surface is likely to be very significant. 


(29) 


(30) 


3.2 Radial Light Rays & the Spacelike/Timelike Conditions 


For radial null paths we have, ds 2 = d6 2 = dcj) 2 = 0, and thus, (1) leads to 


dt R! 

d r n \A + /’ 


( 31 ) 


where the + & — signs are used for outgoing and incoming paths respectively. It is clear from section 
2.5 that the appropriate expression for R! depends on the value of /, so we must consider elliptic and 
hyperbolic cases seperately. The ratio of the null slope (31) to the slope of S(r) (30) is 2 


n 


\R'\ 

a 'V i + /1 


(32) 


Locally, the surface S(t) is timelike if it is steeper than the (radial) null path, and spacelike if it’s less 
steep: 


S(t) timelike: 


o < n < i 


1 Eq (26) follows from R ab = ( T ab - (1/2 )Tg ab ) A 2 = A, - (1/2)T. 

2 We only need the ratio of absolute values, since the slopes of the incoming and outgoing light rays (31) have 
the same magnitude. It is convenient to put the simpler expression in the denominator. 
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S(t) null: 

S(t ) spacelike: 


n = 1 

1 < 7Z < oo 


( 33 ) 


Generically, owing to ( 12 ), R'/\J 1 + / y4 0 even at an extremum, only approaching it near the bang 
or crunch. Therefore, a simultaneous bang, a' = 0, ensures 1Z = 00 and S(t) everywhere. For all 
other considerations, the hyperbolic and elliptic cases need to be considered separately. 

We wish to compare the location of any timelike S(r) region with the horizon loci of the model. 
The apparent horizon (AH) is the locus along a radial light ray where R changes from increasing to 
decreasing. Using equations (31) and (4) this gives 


R'i 

which leads to 



dR 

dr 


= 0 = R 


d t 
dr 


+ Rf = R 


±R' 

\/l + / 


+ & = 


Ty/2 M/R + f 

±\/l + / 


+ 1 


By equations (7) and (9) 

cosh tiah = 2/ + 1, 
COS (rj AH ) = 2/ + 1, 


Rah = 2 M. (35) 

tAH = a+ (sinh rj AH - r] AH )M // 3/2 , (36) 

t A H = a + (rj AH - sin r) AH )M /(—/) 3/2 , (37) 


for hyperbolic and elliptic cases, respectively. 


3.2.1 Hyperbolic Regions 


Substituting (18) into (32), we get 


n 


1 

M' 

W \/i + / 

f 


Mf 

+ 




(38) 


The relative sizes of the ^ functions varies greatly throughout the evolution, therefore we separate the 
model lifetime into three different regions, early, intermediate and late, and we investigate the ratio 
(38) in each of these regions. The early and late time behaviour of the various evolution functions ^ 
were summarised in 52.5. 


• At early times (r/ — > 0), (38) reduces to 


n 


i 'W 1 + / 


M'r? Mfrf r- 2 

TI + ~Fm + ( ~ a) ' /f v 


f 


V V 1 + /’ 


(39) 


assuming a' yf 0, and it is clear that (39) goes to infinity. This means that the surfaces S(r) 
become spacelike at early times. Thus, for any a', the S(t) are spacelike everywhere at early 
enough times. 

• At intermediate times we have to deal with the full expression (38), so only a rough indication 
can be obtained. It makes sense to look for the minimum of 1Z as 77 changes, dlZ/dq = 0. The 
derivatives of ^ 4 , -0 5 and '06 are given in (22)-(24) and plotted in figure 1. 



4- 


V 

2 - 



Figure 1: The shape of d^/d rj, dijj 5 /dr] and d'0 6 /dr/ for the hyperbolic case. 


These 77 functions are multiplied by the same functions of r as in (38), so the exact minimum 7] 
will be different for each worldline. But if for example we assume the r functions are of similar 
magnitude, then r] min fa 1.7. Now, (38) reduces to 

^ \AfM' + BMf' + C(-a')f 5 / 2 
|a'|/ 2 x/T + 7 



where A = i>4{v m i n )- 5 = i’siVmin) and C = 7jj 6 (r] min ) are constants. At r/ mm = 1.7 they are 
A = 0.46, = 0.22, C = 1.45. For S'(r) to be timelike at rj = r] min , 1Z < oo in (40) leads to 



T 

a' 


M' f 1 

M +B J 



(41) 


where T is the scale time of (10). It would appear to be quite easy to satisfy this constraint by 
making |a'| large enough and M' & f small enough. 


At late times (77 —> 00), (38) reduces to 

77 

If a' 7 ^ 0, timelike S(r) requires 


1 

M' 

\ a ' h/l + / 

f 


M f e v . ,, n 

+ ~~p 4" + (~ a )v / 7 


M' = /' = 0, 


(42) 

(43) 


in which case (42) reduces to 77 = \Jf /(l + f) < 1, and this is always true, because / > 0. 
Generally, S(r) is timelike at late times only if M' — f — 0 and a! ^ 0 , otherwise it is spacelike. 


3.2.2 Elliptic Regions 

Using (14), equation (32) evaluates to 

1 

77 = 


Klv/i + / 


M' M f' , - 

+ (—«')>/(—/)^3 


(44) 
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At early times (77 — > 0), (44) reduces to 


TZ = 


?/|\/l + / 


M' 


rj 


Mf f 

(-/) 6 ' 7^40 


+ 


+ 


-a')VFn- 



i + /’ 


(45) 


assuming a' f 0. Clearly (45) diverges, which means S'(r) spacelike at early times for a' f 0, as 
well as a' = 0 . 


• At an intermediate time, which we choose to mean maximum expansion, 77 = tt, equation (44) 
reduces to (see table 2 ) 


1 - 

2M'f + 2Mf\ 


1 

- 2 M' | 2 Mf 

2\R!\ 



M/Vi + / 


WW 1 + / 

f + r- 

MVi + / 


(46) 


where R' is the r derivative of the scale length defined in (10). Applying the timelike condition 
TZ < 1, (46) gives 


\a > 


\R'\ 


VT+7' 

As with the hyperbolic case, this is quite easy to satisfy on any given worldline. 
At late times (77 = 2n — 5, 5 —> 0), if a' ^ 0, equation (44) reduces to 

1 


1Z = 


|a'|\/l + / 
2|27rT' + a' 


M' 4t r Mf 6 tt . f . n—^(- 2\ 

T + 73-T + ( "“ 1 CCl) {—) 


-/) 


|a'|<5 


1 + /’ 


(47) 


(48) 


T' being the r derivative of the scale tine in (10). It is clear that TZ diverges unless 2nf' + a' = 0, 
which makes TZ = 0 and S'(r) a timelike surface. We find that S(r) becomes timelike near the 
big crunch only if the crunch time is constant and the bang time is not. 


Finally, we can summarize these results in table 3. 


4 Explicit Models 

We considered several different time regimes of Lemaltre-Tolman models in section 3, and derived 
the necessary and sufficient conditions for the constant Ricci time surfaces, S'(r), to be timelike or 
spacelike. The main purpose of this section is to investigate numerically the detailed behaviour of these 
surfaces in a range of different explicit LT models, including elliptic and hyperbolic cases, and thereby 
visualize the character of S'(r) in (t, r) plots. We also check if the regions of timelike S'(r) occur inside 
or outside the apparent horizons. The models considered below are only intended to be realistic over 
the finite range of r that is plotted. 
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Early times 

Intermediate times 

Late times 

Hyperbolic 

Always 

spacelike 

S(t) timelil 

and a' ^ 0, 

:e only if 

TiKA 

otherwise spacelike 

S(t) timelike if 

M' = f = 0 

and a' yf 0, otherwise 
spacelike 

Elliptic 

Always 

spacelike 

S(t) timelike only if 

, „ \R’\ 

Vi + / 

and a' ^ 0, otherwise spacelike 

S(t ) timelike if 

2 v tT ' + a' = 0 

and a' y 4 0, otherwise 
spacelike 


Table 3: Summary of the conditions for S(t) to be timelike or spacelike at early, middle and 
late times. It is understood that S(r) is null on the boundary between timelike and spacelike 
regions. 


4.1 Numerical Calculations 

We have developed a numerical code that takes the definitions of the LT arbitrary functions and 
calculates dt/dr\ n , dt/dr| T and 1Z, as detailed in section 3. The resulting data is plotted against r and 
t as a pair of surfaces of different colours, whose heights above the (t,r) plane are the slopes df/dr| n 
and dt/dr\ T . In order to check that our model choices were reasonable and free of irregularities, the 
code also produces plots of M(r), f(r), a(r ) and T(r), as well as the functions that indicate a shell 
crossing if they have the wrong sign, M'(r), f'(r), a'(r ) and 2nT'(r) + a'{r). Out of these last 4 
functions, the first 3 are relevant to hyperbolic regions, and functions 1, 3 and 4 are relevant to elliptic 
models (see table 1). 

The actual code must allow for a few practical numerical issues. 

• The slope df/drj n diverges near the bang and the crunch, and this must be chopped by not 
plotting certain q values. 

• In elliptic regions, if the big crunch occurs within the time range plotted, the two surfaces must 
be terminated just before the crunch. 

• In calculating the locus of the AH, it is easiest to use (36) & (37), since we are already using q 
in the evolution calcuations. In fact there are two AHs in elliptic regions. Although each AH is a 
locus in the (t,r) plane, for clarity, the AH curves are plotted on the null slope surfaces, i.e. at 
“height” df/dr| n on the plot, otherwise they would be hidden. 

• As noted above, the origin and spatial extrema in R generate undefined zero-over-zero values 
unless their limits are evaluated properly. Therefore, the various arbitrary funtion combinations 
appearing in equations (38) and (44), such as Mf'/f 2 /y/1 + f/a' and f/{l + /), must be 
re-evaluated analytically for each choice of M, f and a, and written into the code explicitly. 

The general algorithm is as follows. 

1. Define the ranges of r and q (typically r changes from 0 to 0.9, while q extends from 0 to 7.2). 
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2. Define the arbitrary functions M, f, a and their derivatives M', /', a! , as well as other ratios 
like M/f and M// 3 / 2 . 

3. Rescale 77 to be fj given by (25) in order to avoid vastly different timescales in the same plot, 
and handle near-parabolic regions without trouble. 

4. In terms of fj, calculate the ^ for elliptic regions from (15), (16) and (17), and for hyperbolic 
regions from (19), (20) and (21). Redefine these functions as series approximations for use when 
-q approaches zero or 2n. 

5. Depending on the value of fj, calculate R, t, dt/dr\ n , dt/dr | T and 1Z for the whole range, using 
the appropriate elliptic or hyperbolic expressions. 

6 . Plot the two slopes as intersecting 2-d surfaces in a 3-d plot . 3 

7. Evaluate the AH locus, and draw it as a black line on the dt/dr \ n surface. 

8 . Plot the LT arbitrary functions and the no shell crossing condition functions. 

4.2 Hyperbolic Models 

4.2.1 Model 1 

The arbitrary functions 


M = M 0 (r 3 + Mir 4 + M 2 r 5 ), (49) 

/ = -k(r 2 + /yr 3 + / 2 r 4 ), (50) 

a = a 0 + ayr + a 2 r 2 + a 3 r 3 , (51) 

describe an inhomogeneous model with non-zero density everywhere, so it could be part of a “cosmol¬ 
ogy”. The first set of coefficients is M 0 = 1, Mi = —3, M 2 = 2.4, k = —1, f\ = —2.67, / 2 = 2, 
a 0 = 5, Oi = —200, a 2 = 0.5 and a 3 = 0.9. These functions and coefficient values were specifically 
selected so that both M and / have zero gradient at the same r value, but never have a negative 
gradient, as shown in figure 2a. Therefore we have a worldine (r value) which satisfies the condition 
(«), M' — f = 0, which should make S(r) timelike at late times, and the neighbouring worldlines 
are very close to that condition. Figure 2a plots the arbitrary functions M, f and a. Scale time T is 
not necessarily increasing with r, but has a zero gradient at the same point as functions M and /. 
Figure 2b shows the conditions for no shell crossings are satisfied, because M' > 0, f > 0 and a! < 0 
clearly hold. The fourth condition 2nT + a > 0 is not relevant for hyperbolic models. 

Figure 3 illustrates the relation between the {constant r} surfaces S(r) (red) and the radial null 
surfaces (blue) for each worldline and a sufficiently large range of the evolution. The black line shows 
the apparent horizon, which allows us to see whether the timelike S(r) regions occur inside (at smaller 
R) or outside (at larger R) the AH. Our hyperbolic model is always expanding, so smaller R occurs to 

the past of the black line. For the worldlines (0.3 < r < 0.7), the S(t) are spacelike at early times, 

they become timelike at intermediate times and remain so for the rest of their plotted evolution. This 
confirms that S(r — > oo) is timelike on r = 0.5 for which M' — f — 0. For the remaining worldlines, 
(0 < r < 0.3) and (0.3 < r < 0.7), the S(t) are spacelike at early and late times, but change to 
timelike for at least part of the evolution. We noticed that the size (duration) of the intermediate 

3 Several other plots, such as 1Z against t and R, are plotted but not used here. 
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timelike region, can be increased by increasing (—a 7 ). Conversely, decreasing (—a 7 ) makes it shrink 
towards the bang. This model confirms that S(t ) is spacelike everywhere for early times, but at late 
times S(t) can be timelike only if M' and f are zero while a' is not. It also illustrates that timelike 
(r) at intermediate times is not hard to achieve. The black line indicates that the timelike region is 
located outside the AH, i.e. at larger R. 



Figure 2: Arbitrary Functions for Hyperbolic model (1) using the first set of coefficients. The model 
is free of shell crossings if the M' and f graphs are never negative, and the a' graph is never positive. 
See table 1. 



Figure 3: Hyperbolic model (1) using the first set of coefficients: The slopes of the S(t) and radial 
null surfaces, dt/dr | T and dt/dr\n, shown in red and blue respectively, are plotted against t and r. 
Where the red surface is above the blue one, S(r) is timelike. The black line represents the apparent 
horizon, and the big bang is where the blue surface rises at the back. 
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Now, by changing the coefficients of the arbitrary functions in equations (49)-(51), we can produce 
a model for which S(r) is spacelike everywhere. Taking M 0 = 200, M\ = 0.5, M 2 = 0.5, k = —1, 
fi = 0.6, / 2 = 0.6, ao = 5, ai = —1, a 2 = —0.5 and a 3 = —0.9, we get the LT functions shown 
in figure 4a, and figure 4b shows there are no shell crossings. The blue dt/dr\ n surface in figure 5 is 
everywhere above the red dt/dr\r surface, showing that S(r) is always spacelike for each world line. 
Th is is because M' > 0 and /' > 0 and a' is not so large. The bold black line visualizes the apparent 
horizon. 

By using a third set of arbitrary functions and varying just one of them, we can investigate the 
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(a) The behavior of M, /, a and T (b) Conditions for no shell crossings 

Figure 4: Arbitrary Functions for Hyperbolic model (1) using the second set of coefficients. 



r Go 


Figure 5: Hyperbolic model (1) using the second set of coefficients: The slopes of the S'(r) and radial 
null surfaces are plotted, showing S'(r) is spacelike everywhere. 
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nature of the intermediate time. Choosing M 0 = 1, Mi = —2, M 2 = 1.4, k = —1, fi = —1.65, 
f 2 = 1 , a 0 = 0 , a\ = — 200 , a 2 = 0.1 and 03 = 0 . 1 , the behavior of the functions and the conditions 
for no shell crossings are represented in Figures 6 a and 6 b, respectively. The coefficient a\ will be 
varied. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


E-h 





r 


(a) The behavior of M, /, a and T 


(b) Conditions for no shell crossings 


Figure 6 : Arbitrary Functions for Hyperbolic model (1) using the third set of coefficients and a\ = —200. 



Figure 7: Hyperbolic model (1) using the third set of coefficients and = —200: The slopes of the 
S(t) and null surfaces are shown, and here S(r ) is timelike only at intermediate times. 

With a! = —200, which creates a moderately large value of (—a'), the S(t ) are spacelike for 
all r, at both early and late times, but they still become timelike at intermediate times, figure 7. If 
(— a') is made very large, a\ = — 2000 , the S(t) are of course spacelike at the earliest times, but 
become timelike just after the AH, and remain timelike up to very late times on all worldlines, figure 
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8 a. However, since M' > 0 and /' > 0, the S(t) must eventually become spacelike again. If (—a') is 
made small, a± = —0.2, the S(t) now become spacelike everywhere at all times, figure 8 c. In other 
words, by decreasing (—a 7 ), the extent of the timelike region can be decreased and made to vanish for 
small enough (—a'). Moreover, the timelike region shrinks towards the origin as (—a') is decreased. 



Figure 8 : Hyperbolic model (1) using the third set of coefficients and 3 different a\ values: This is a 
view looking straight “down", so the “heights” dt/dr are not seen. The sizes of the red regions show 
the sizes of the regions of timelike S(t). The black line is the apparent horizon. 


4.2.2 Model 2 

The arbitrary functions are 


M = M 0 exp(l + Mir), (52) 

f = —k exp(l + fir), (53) 

a = ao exp(l + air), (54) 

where M 0 = 2, M x = 0.01, k = —1, f\ = 0.1, a 0 = —0.01 and a\ = 5, and they’re plotted in figure 
9a. This is a hyperbolic model with an unusual topology, because the “origin” is asymptotically far 
away at r = — oo. However, figure 9b shows it is free of shell crossings. 

At early times, the S(r) are spacelike everywhere, but become timelike at intermediate times for 
a specific region. As before, the size of the timelike region depends strongly on (—a 7 ), and it occurs 
outside the AH. At late times, the S(r) eventually become spacelike again. 

4.3 Elliptic Models 
4.3.1 Model 1 

The following arbitrary functions 


M = M 0 (r 3 + M x r 4 + M 2 r 5 ), (55) 

f =-k(r 2 + fir 3 + f 2 r 4 ), (56) 

a = a 0 + air + a 2 r 2 + a 3 r 3 , (57) 

are first considered with M 0 = 1, Mi = 0, M 2 = 2.4, k — 1, fi — —0.67, f 2 = 0, ao = 5, ai = —5, 

a 2 = 0.5 and a 3 = 0.9, as plotted figure 11a. This function set represents an inhomogeneous elliptic 
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model, which could be part of a cosmology. Although these functions would not work globally, they 
are entirely valid for the range of r calculated. Figure lib shows no shell crossings occur; for elliptic 
models the /' plot is not relevant. 

Figure 12 illustrates the relationship between the S(r) and the null surfaces, whose slopes are 
represented by the heights of the red and blue surfaces respectively, in the 3-d plot. The bold black 
lines show the loci of the apparent horizons. At early and late times, the S(t ) surfaces are spacelike 
everywhere, but the S(r ) become timelike at intermediate times, for a range of worldlines (0 < r < 0.3). 
The extent of the timelike region once again depends on the gradient of the bang surface, and as a' 
decreases the region near the origin disappears last. The region of timelike S(r) lies outside the two 



(a) The behavior of M, /, a and T 
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(b) Conditions for no shell crossings 

Figure 9: Arbitrary Functions for Hyperbolic model (2). 



Figure 10: Hyperbolic model (2): The slopes of the S(t) and null surfaces are indicated by the heights 
of the red and blue surfaces in this plot, and S(r) is timelike where the red surface is higher. 
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AHs (larger R), though very close to them for quite a range of r. 

To investigate how the size of the timelike region is affected by the crunch time, we next keep the 
scale time T unchanged, but adjust a to make the crunch time constant. Thus M(r) and f(r) are as 



(a) The behavior of M, /, a and T 




(b) Conditions for no shell crossings 


Figure 11: Arbitrary Functions for Elliptic model (1) with the first coefficient set. The model is free of 
shell crossings if the M' and 2tcT' + a' graphs are never negative, and the a' graph is never positive. 
See table 1. 


100 , 



Figure 12: Elliptic model (1) with the first coefficient set: The relationship between the 5(r ) (red) 
and null (blue) surfaces, where their heights represent the slopes of these surfaces. The plotted time 
range only includes the crunch at smaller r values. The edge of the blue surface is irregular due to 
the clipping of divergence near the bang and crunch. The black lines represent the past and future 
apparent horizons. 
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in (55)-(56), but a becomes 


a = -27rT = -27T-—, 
J 6/z 


( 58 ) 


The coefficients are changed to M 0 = 1, Mi = 0.1, M 2 = 0.1, k = 1, f± = —0.67 and / 2 = 0. See 
figures 13a and 13b. 



0i-1-1- f-t -1-1-1-1- 0 



0 I-1-1-1-'-'- 1 - 1 - 1 - CM- il -1-'-'- 1 - 1 - 1 - 1 - 1 - 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


r r 

(a) The behavior of M, /, a and T (b) Conditions for no shell crossings 

Figure 13: Arbitrary Functions for Elliptic model (1) with the second coefficient set. 



Figure 14: Elliptic model (1) with the second coefficient set: The relationship between the S(t) and 
null surfaces. The crunch is not included in the plotted range for the larger r worldlines. 

Figure 14 shows the S(t) surfaces are initially spacelike everywhere, but at intermediate and late 
times they become timelike, owing to the constant crunch time. 
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This situation can be changed if a is made a different multiple of lifetime T. For instance, if 
a = —0.9(27rT), the S(t) are spacelike at early times, and stay timelike through intermediate up to 
quite late times. If a = —0.2(2 tcT), the S(t) become timelike at intermediate times and return to 
spacelike quite soon thereafter, as shown in figure 15. In this case, we find two discrete timelike regions; 
one close to the past AH, the other is close to future AH and the origin. If a = —0.01(27rT), the S(r) 
become spacelike everywhere. As always, the timelike region or regions are always outside the apparent 
horizons (the black lines), though often very close to them. 



Figure 15: Elliptic model (1) with the second coefficient set and a = —0.2(2 ttT): Comparison of the 
S(t) and null surfaces, showing the S(r) are timelike only at intermediate times. The plotted evolution 
range does not reach the crunch at larger r values. 


4.3.2 Model 2 

For this model the arbitrary functions are 

M = M 0 + Mi exp(r 2 /M 2 ), (59) 

/ = -A;exp(r 2 // 2 ), (60) 

a = a 2 r 2 , (61) 

where M 0 = 2, Mi = 1 , M 2 = 0.1, k — 1, f 2 — 0.1 and a 2 = —1; see figures 16a and 16b. 

Th is set of functions represents a non-vacuum black hole; it has /(0) = —1 and /'(0) = 0 = 

M'(0) = a'( 0 ), as required to get the full manifold with two asymptotic regions joined by a “neck” 

[18]. There are two asymptotic regions, and matter flows out of the past singularity (or bang) and into 
the future singularity (or crunch). Figure 17 shows the S(t ) surfaces become timelike at intermediate 
times, outside the black hole, while at early and late times the S(r) are spacelike everywhere. A 
vacuum black hole is obtained by putting M' = 0, but note that it is not appropriate for testing Ellis’s 
idea because there are no unique timelike Ricci eigenlines. 
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4.4 Discussion 


We have found in all cases that the {r= constant} surfaces S(t) are spacelike everywhere at early 
times. 

For the hyperbolic models examined, timelike regions seem common at intermediate times unless 
a' is very small. However the S(t c ) are nearly always spacelike in the late time limit, except for the 
rather special case that M' — f — 0 on certain worldlines. For worldlines that are very close to this 
condition, the timelike character of the S(r) can remain until for a long period. Now M' = 0 ensures 



(a) The behavior of M, /, a and f (b) The no shell crossing conditions 

Figure 16: Arbitrary Functions for Elliptic model (2); a non-vacuum black hole. 



Figure 17: Elliptic model (2): The character of the S(r) and null surfaces for a non-vacuum black 
hole. 
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zero density locally, but if it’s only at one r value, the Ricci eigenline still exists as a limit. 

For elliptic cases, intermediate times could easily have timelike S(t) regions, but to get timelike 
S(t) at late times (the crunch), the crunch time must be a constant, and the bang time not. It is 
notable that with black hole and black-hole-like models, the regions of timelike S(r) do not occur along 
worldlines near to the “neck” but they do occur further out. 

In the models looked at, the timelike regions always occur outside the horizon, but often very close 

to it. 

Invariably, the size of a timelike region was increased by increasing {—a') and vice versa, vanishing 
completely for samll enough (—o'). Additionally, as {—a') decreases, the timelike region shrinks towards 
the origin. 


5 Conclusion 

The true nature of the past, present and future is a mystery, and how the flow of time can best be 
represented diagrammatically remains problematic, so there has been a great effort by many authors 
attempting to envisage the flow of time and give a clear picture of our universe [8, 1, 33]. According to 
the block universe idea [14], the universe consists of a unique family of surfaces Sir), and the present is 
merely one surface of constant r. Conditions at the present time can be evaluated by integration from 
known conditions at any earlier time, possibly the beginning, up to the present time. The evolving block 
universe (EBU) and the crystallising block universe (CBU) take the second law of thermodynamics and 
quantum effects seriously, and replace the time-reversible, deterministic model with one that recognises 
the past is fixed while the future is highly mutable [13, 15]. 

Our main purpose was to investigate the viability of a particular proposition of Ellis [14], that the 
correct time r to use is the proper time along the timelike eignevectors of the Ricci tensor (“Ricci 
time”). It was suggested that if the resulting constant r surfaces “sensed” strong gravitational fields, 
for example by becoming timelike, this would indicate the right kind of behaviour. 

The Lemaltre-Tolman (LT) metric [22, 27, 17] was used to generate a variety of models with 
different features; some like inhomogeneous cosmologies, and some like matter-filled black holes. The 
slopes of the constant r surfaces S(t) were compared with the slopes of the radial null surfaces in 
order to determine if they were spacelike, null or timelike. 

Analytic methods were used to find conditions, on the LT arbitrary functions, M(r), f(r) & a(r), 
for the Sir) to become timelike in two limits: early times near the bang, and late times near the 
crunch for elliptic models or r —» oo for hyperbolic models. Intermediate time calculations, though not 
definitive, indicated it was rather easy to make S(r) timelike. Numerical methods were then used to 
calculate and display the full evolution of this relationship for a variety of explicit LT models. 

In summary, the S(r) surfaces are spacelike everywhere at early times for models that are free of 
irregularities; it is already known that the bang surface is spacelike. In the late time limit, the S(r) 
can be timelike, though only under quite strong conditions, and these conditions are totally different 
in hyperbolic and elliptic models. In the hyperbolic case, both the geometry/energy function f(r) and 
the interior mass function M(r) must be (locally) constant, M' = 0 = /'. This means there is vacuum 
p = 0 on these worldlines. In the elliptic case, the crunch time must be constant (as it is in FLRW 
models), while the bang time is not constant. 

We found that timelike S(t) regions are common at intermediate times. The “area” of the timelike 
S{t) regions (in the (t,r) plane) strongly depends on (—a'). The range and duration of the timelike 
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S(t) regions is increased by making (—a') bigger, and if (—a') is sufficiently small the S(t) surfaces 
will be everywhere spacelike. Also, as (—o') is decreased, the timelike regions sometimes tend to shrink 
towards the origin (if there is one). 

We note that the condition, a' = 0, i.e. a constant bang time, which guarantees all constant Ricci 
time surfaces are spacelike throughout a model, is also the condition which ensures there are no decaying 
modes in LT models [38]. The growing and decaying modes have been invariantly characterised in [41]; 
see [35] for the definitions. The asymptotics of LT models, represented by covariant scalars along radial 
geodesics, including those with a' = 0, were investigated in [40]. Of interest here is the relation to 
gravitational entropy, which was also investigated for LT models in [43]. They considered the entropy 
proposals of [10, 25] and a modification of the latter in [42], They found that entropy grows if there 
are anti-correlations between Hubble and density fluctuations. This requires that the decaying modes 
are sufficiently suppressed, and near the bang they have to be absent for entropy to grow. They 
also suggest that regions of decreasing entropy and dominant decaying modes may be associated with 
instability. A related result [5] is that for generic inhomogeneous models, the inhomogeneity initally 
decreases and then increases, and that entropy decreases while inhomogeneity decreases. 

The timelike S(t) regions always appeared outside the apparent horizons, i.e. at larger areal 
radius, though they were often right next to them. 

For one specific elliptic model we found there were two discrete timelike S(r) regions on certain 
worldlines, separated by a significant time. One of them was near the past apparent horizon (AH) and 
the other was very close to the future AH. 

In the black-hole-like models, we found that the timelike regions not only stayed outside the AHs, 
but also did not occur on worldlines at or near the “neck”. 

Overall, we found that the Ricci time has the opposite of the suggested behaviour, and if its 
{constant r} surfaces S(r) become timelike, this is in regions outside the apparent horizons, where 
the gravitational fields are not so strong. Therefore this time may not be appropriate for the EBU. 
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